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Evidence that the Y(4660) is a / (980)i/>' bound state 
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Abstract 

We demonstrate that the experimental information currently available for the Y(4660) is consistent with its being a /o(98Q)V>' 
molecule. Possible experimental tests of our hypothesis are presented. 
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1. In recent years a large number of new hidden and 
open charm states were discovered experimentally above 
the first inelastic thresholds. For most of them the masses 
have in common that they deviate significantly from the 
predictions of the quark model [l| — which on the other 
hand works very well below the thresholds — and are 
positioned very close to a s-channel threshold For 
recent reviews see, e.g., Ref. 0]. 

The proximity of the thresholds initiated a lot of the- 
oretical studies, in order to reveal a possible molecular 
structure of those states. E.g. in Refs. [4| the £> s (2317), 
located just below the i^D-threshold, was studied within 
the molecular model and the V(3872), located right at 
the DD* threshold, was investigated in Refs. [5J. How- 
ever, so far no consensus exists on the true nature of 
those states and, e.g., four-quark interpretations [6j and 
even conventional qq states are still under discussion 0] . 

In Refs. it was argued that there is a way to model 
independently identify hadronic molecules in the spec- 
trum. The method is based on a time-honored analysis 
by Weinberg [§] and applies, if the pole is very close to 
the threshold of the constituent particles that form the 
bound state in an s-wave. The method relates the ef- 
fective coupling constant of the bound state to its con- 
stituents, g, directly to the molecular admixture of the 
state. Especially, one may write for a pure molecule 
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where mi and ttt-2 denote the masses of the constituents 
and e the binding energy. 

In this paper we discuss the nature of the Y"(4660) 
as a candidate for a fo(980)ip' bound state. So far 
the y(4660) was seen only in the 7r + 7r _ -0' invariant 
mass distribution in e + e~ — > 7/s_R7r + 7r _ i// with a mass 
of 4664 ± 11 ± 5 MeV and a width of 48 ± 15 ± 
3 MeV [l(|. Such a structure was neither observed 
in e + e~ — > 7/s , fi7r + 7r _ J/ip [11], nor in the exclusive 

nor 
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e+e- -> DD,DD*,D*D*,DDtt cross sections 
in the process e+e~ J / 'ip Z)M £>W pi]. These facts 
would severely challenge any attempt explaining the state 
as a canonical cc charmonium, e.g. 5 3 <Si as in Ref. 
14j . The Y74660) is suggested to be a baryonium state 
in Ref. [15j|. The difficulties in the canonical char- 
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FIG. 1: Diagram illustrating the most prominent decay chan- 
nels of the F(4660) in the molecular model. The solid black 
dot denotes the Yfoip' vertex, whose strength parameter g is 
fixed within the molecular model through Eq. |T}. The de- 
cay vertices for /o, shown as open circle, are fixed from other 
data. 



monium interpretation could be explained naturally, if 
the y(4660) were a /o(980)^' bound state, because it 
would decay dominantly via the decay of the /o(980), 
i.e. F(4660) -> V7o(980) -> ^'[tttt] and F(4660) -> 
V>'/o(980) -> i/j'[KK]. The nominal threshold of the 
/o(98O)--0' system is at 4666 ± 10 MeV if we take the 
PDG value of the /o(980) mass 0. Note that the inter- 
action between a heavy quarkonium and a light hadron 
is expected to be dominated by the QCD van der Waals 
interaction which is attractive (l7| . 

In the 7T7T invariant mass a clear /o(980) peak is vis- 
ible. Clearly, if the F(4660) were a conventional foifj' 
bound state, it could not decay to this channel. This is 
possible only because of the finite width of the /o(980), 
mainly due to its decay to the 7T7t channel. We will per- 
form our analysis based on the following reasoning: if the 
/o were a stable particle, also the F(4660) were stable. 
Then we could calculate the effective coupling constant 
of Y to fo using Eq. (fl]). The central assumption is 
that this coupling constant does not change as the two 
pion channel is switched on — a similar ansatz lead to 
a successful phenomenology for the /o(980) treated as a 
KK molecule 18]. Under this assumption we can pre- 
dict not only the prominent component of the width of 
the F(4660) but also its spectral shape using the mass 
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and the over all normalization as the only input. As addi- 
tional non-trivial result we can also predict the strength 
and the shape of the decay to KK and 77. The possi- 
ble decay chains are illustrated in Fig. [TJ As we will see, 
the resulting spectra are in excellent agreement with the 
existing data. We interpret this as a strong support in 
favor of a molecular interpretation for the F(4660). 

2. The y(4660) was observed as a structure in 
the iP'tttt invariant mass distribution. In addition, in 
Ref. [Toj ] . also the irir invariant mass distribution is pre- 
sented, which was found after applying an appropriate 
cut to the ip'iriT invariant mass. Both distributions can 
be derived from the same differential rate, under the as- 
sumption that the pion pair stems from the decay of a 
/ . One gets 
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also [22j). The input parameters are taken from the fits 
provided in Ref. [19(. To be concrete, for all the curves 
shown below we used the central values of the various pa- 
rameters of fit K2 shown in table 4 of that reference, thus 
we used m/ = 0.9862 GeV, gf x+K- = 3.87 GeV, and 
S/ott+tt- = —2.03 GeV, while the couplings for the neu- 
tral channels were fixed using the isospin relations. We 
checked that the other parameter sets give very similar 
results to the ones discussed in detail below. 

The only missing piece is the physical Y propagator. 
In Eq. ([3]) an explicit expression is given for the par- 
tial width to t//7t + 7t~ . However, in order to derive the 
physical propagator, all relevant decay channels need to 
be included. Within our model we assume that the Y 
decays predominantly through the fo, thus a consistent 
treatment requires that all decay channels of the fo also 
contribute to the width of the Y and thus to the prop- 
agator. It is straightforward to extend Eq. (J3]) also to 



where Gy{M) denotes the physical F-propagator, to be 
specified below and iV is a normalization constant, which 
contains, besides the electron-photon vertex and the pho- 
ton propagator, both indeed constant to very high accu- 
racy in the range of M of relevance here, also the detector 
acceptance. The decay Y — > ip'ir + Tr~ is described by 
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Here p denotes the c.m.s momentum of the tp' for given 
values of the irir invariant mass, m^, and the ip'mr in- 
variant mass, M, and M^/ W7r = m^r + M^,/. The effective 
coupling constant g is not a free parameter, but can be 
determined for any given value of the mass of the Y from 
Eq. |T]). The quantity w \m-nir) denotes the tt + it~ 
fraction of the fo spectral function. It is normalized ac- 
cording to 
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any additional free parameters — also for these channels 
we use the results of the fit to the KLOE data presented 
in Ref. [lj|. Thus we get for the total width of the Y 
at a given value of M, under the assumption that it is 
saturated by the tp'fo decay 
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where s^ r 



(m a + nib) 2 ■ In order to get a Y propagator 
with the correct analytical properties we need to continue 
the contribution from the KKip' channel also to below 
its threshold. For this we use a dispersion integral, which 
gives us an expression for the Y self-energy, Hy(M), for 
arbitrary values of M 



7T Jm* s- M z - ie 
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With this normalization Eq. ([3]) goes to the standard 
expression for a two particle decay in the stable particle 
limit for the Jq. 

A high quality data set for the /o was collected recently 
by KLOE [19( based on the reaction <j> — > jirir. The data 
was analysed using the so-called kaon loop model [2(| 
and provided parameters for the /o with very little un- 
certainty. To be concrete, we use 
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where the ft^(m) = ^f {m) - Re(H^(m /o )) denote the 
renormalized self-energies of the fo with respect to the 
channel ab. Analytic expressions are given in Ref. [2(| (see 
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2m^ denotes the lowest physical 



threshold of relevance here. Note, this treatment is com- 
pletely consistent to what was done for the fo — one way 
to derive the self-energies n^(m) given above is through 
a dispersion integral with the two-body phase space as 
input. With the self energies at hand we may now give 
the expression for the physical propagator of the F(4660) 



G Y (M) = 



M 2 - M Y +U Y (M) 



(7) 



where, as above, we defined Ily(M) = Hy(M) — 
Re(ny(Mr)). 

3. Our model has only 2 free parameters, namely, the 
mass of the Y, My, and the normalization constant, N, 
introduced in Eq. @. We now proceed as follows: the 



2 




4.6 4.7 4.! 

M(v|/tiV ) (GeV) 



4.9 



FIG. 2: Comparison of the line shape in the ip'iv + Tv~ invariant 
mass distribution derived from the molecular model with the 
data. The solid line shows the result of our best fit, while the 
shaded area shows the uncertainty that emerges from the fit. 
The dashed line shows the best fit result, if also the effective 
coupling is part of the fit, as discussed in the final section. 



count rate, R, in the ip'^ir invariant mass spectrum is 
given by 



R(M)-- 



dm 2 
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(8) 



using the expression for the differential rate W of Eq. 
The range of parameters allowed by the data is then 
determined from a fit to the experimental data for the 
iP'tt + tt^ spectrum — we remind the reader that the cou- 
pling constant g is fixed from Eq. |1| for any given value 
of M y . We find 



N = 10 ± 2 GeV 3 , M Y = 4665t? MeV 



(9) 



This 



range of mass values corresponds to a range of 
g = 11... 14 GeV for the effective coupling constant. 
The best fit is shown as the solid line in Fig. [5] — the un- 
certainty that emerges from the fit is reflected by the grey 
band. For the best fit x 2 /d.o.f = 0.5. The first observa- 
tion is that the resulting invariant mass distribution visi- 
bly deviates from a standard Breit-Wigner shape. This is 
a direct consequence of our starting assumption, namely 
that the F(4660) is predominantly composed of a /o and 
a ip' , not only since the mass of the Y is very close to the 
nominal foip' threshold, but also because of the proxim- 
ity of the kaon channels, which are very important for 
the structure of the /o- 

Since now all parameters of our model are fixed we can 
predict other channels. Recently high quality data was 
measured for fo — > 77 2l| . Within our approach we thus 
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FIG. 3: Predictions for the irn and KK invariant mass dis- 
tributions. 



predict that the shape in the ^'77 invariant mass is iden- 
tical to that measured for , however, down scaled 
by the relevant ratio of branching ratios. In addition we 
can also predict the signals for the spectra of 7T7t and KK 
as they emerge from the decay of the Y. The predicted 
rate, R a b(ni>ab), again follows directly from Eq. ([2]) 
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where the limits of integration are chosen identical to 
the cuts used to get the experimental rate [lOJ. The 
corresponding results are shown in Fig. ©. If a signal 
of the given shape and strength were found in the iji'KK 
invariant mass distribution, it would provide a strong 
support for the assumed prominent role of foip' for the 
structure of the Y. 

4. To summarize, we calculated the invariant mass 
spectrum for i//7r + 7r - as well as the corresponding 7T7t and 
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KK spectra in the mass range of the F(4660) under the 
assumption that the E(4660) is a /o(98O)-0' bound state 
and /o(980)'i/' / is its only decay channel of relevance. A 
very good description of both spectra, where data exist, 
was achieved. Especially, we find a visible deviation from 
a Breit-Wigner shape for the i/j'tvtv spectrum, consistent 
with the data, although data with better statistics would 
be very welcome to strengthen this point. 

In principle our method also allows one to estimate 
the possible contribution of channels other than fy-ip' to 
the width of the E(4660). One just needs to add the 
term iMyT^d to the denominator of Gy(M), defined in 
Eq. (JT]) and repeat the fitting procedure. With this we 
get r a dd = (30 ± 30) MeV — thus, before better data 
is available no reliable bound for the possible additional 
width can be deduced. However, it is important to ob- 
serve that the value extracted from the current data is 
consistent with zero within the uncertainty. 

We also checked what happens, if we do not fix the 
effective coupling g according to Eq. (fTJ), but allow it to 
float as well. Then we find My = (4672 ± 9) MeV and 
g = (13 ± 2) GeV with x 2 /d.o.f = 0.4. The result of the 
best fit is shown as the dashed line in Figs. [2] and [3l Thus 
the effective coupling constant extracted from this three 
parameter agrees to that found before, only that the three 
parameter fit prefers a larger value of the mass, partially 
inconsistent with a molecular picture — for masses at 
the higher end the mass of the Y is even larger than 
Mjji' +rrif . Also here we need to conclude that more data 
are needed to draw a more firm conclusion — especially 
we showed that the KK invariant mass distribution is 
quite sensitive to the mass of the E(4660) — c.f. lower 
panel of Fig. [31 It is in any case important to stress that 
the fit calls for a large coupling constant for Y — > foip', 
which is naturally explained by the assumption that the 
E(4660) is generated dynamically in the foip' channel. 

We take our results as a strong evidence for a molecular 
interpretation of the E(4660). We stress that a measure- 
ment of the tp 1 KK channels as well as an improvement of 
the resolution of the existing data would allow for non- 
trivial tests of our hypothesis. 
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